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Abstract 

We study theoretically the quantum optical properties of hybrid molecules composed of an 
individual quantum dot and a metallic nanoparticle. We calculate the resonance fluorescence of 
this hybrid system. Its incoherent part, the one arising from nonlinear quantum processes, results 
to be enhanced by more than two orders of magnitude as compared to that in the absence of 
the metallic nanoparticle. Scattering spectra at different excitation powers and nonperturbative 
calculations of intensity-field correlation functions show that this system can act as a nonlinear 
ultra-compact two-photon switch for incident photons, where the presence (or absence) of a single 
incident photon field is sufficient to allow (or prevent) the scattering of subsequent photons. We 
also find that a small frequency shift of the incident light field may cause changes in the intensity 
field correlation function of orders of magnitude. 

PACS numbers: 42.50.Ct, 32.70.Jz, 42.55. Sa, 78.67.Hc 
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I. INTRODUCTION 



Control over the interaction between single photons and individual optical emitters de- 
serves great importance in quantum science and quantum engineering^. Recently, substan- 
tial advances towards the realization of solid state quantum optical devices have been made 
coupling single quantum dots (QDs) to high-finesse optical cavities^^^. An inherent limi- 
tation of cavity quantum electrodynamics (QED) is that the size of the cavity is at least 
half wavelength and practically much more than that owing to the presence of mirrors or 
of a surrounding photonic crystal. Unlike optical microcavities, metallic nanoparticles and 
metallic nanostructures are able to focus electromagnetic waves to spots much smaller than 
a wavelength. In this way it is possible to increase the local density of the electromag- 
netic modes as in microcavities but with ultra-compact structures. The ability of metallic 
nanoparticles and nanostructures to control the radiative decay rate of emitters placed in 
their near field has been widely demonstrated^"^. This ability stems from the existence of 
collective, wave-like motions of free electrons on a metal surface termed surface plasmons 
(SP)^. Moreover the plasmon excitation covers a broad bandwidth and requires no spe- 
cial tuning to achieve resonance. An outstanding demonstration of the cavity-like behavior 
of metallic nanoparticles is the recent realization of a nanolaser based on surface plasmon 
amplification by stimulated emission of radiation (spaser)^. 

Optical nonlinearities enable photon-photon interaction and lie at the heart of several 
proposals for quantum information processing-*^'^, quantum nondemolition measurements 
of photons^iii, and single-photon switching^^ and transistorsi^. Recently, the nonlinear 
optical response of a semiconductor QD coupled with a metallic nanoparticle (MNP) has 
been theoretically investigatedi^^i^ by exploiting a semiclassical approach where the quan- 
tum emitter is treated quantum mechanically and the light field classically. Although this 
approach provides useful information on the absorption and elastic scattering of this sys- 
tem, the nonlinear optical properties of individual QDs display important quantum optical 
effect^^~—. Moreover, in order to investigate the possible use of these hybrid artificial 
molecules as quantum devices for the control of individual light quanta, a full quantum 
mechanical description is required. Recently the efficient coupling between an individual 
optical emitter and propagating SPs confined to a conducting nanowire has been studied 
both theoretically^^ and experimentally^^. Non-classical photon correlations between the 
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emission from the quantum dot and the ends of the nanowire demonstrate that the latter 
stems from the generation of single, quantized plasmons. The potential of this system as an 
ultracompact single-photon transistor has been theoretically demonstrated^^. 

Here we investigate the quantum optical properties of a QD-MNP hybrid artificial 
molecule (see Fig. la). We consider a spherical QD interacting with a spherical MNP of ra- 
dius Tm, separated by a distance R (as shown in Fig. la). Throughout the paper we consider 
the exciting input field polarized along the molecule axis (see Fig. la). The quantum dot 
is modeled as a two level system, that is a good approximation when investigating optical 
processes at frequencies resonant with the lowest energy excitonic transition^^. Calculations 
are carried out nonperturbatively, i.e. the QD-MNP interaction as well as quantum fiuctu- 
ations are treated at all orders. We start analyzing the scattering properties of the system 
as a function of incident power. We also investigate the modification of the QD resonance 
fluorescence induced by the presence of the metallic nanoparticle. Resonance fluorescence 
of individual quantum emitters provides an interesting manifestation of the quantum the- 
ory of light^"— . The spectral and quantum statistical properties of the fluorescent light 
radiated by atom emitters have been the subject of quantum optical measurements for 
many years. More recently resonant fluorescence of individual QD has been experimentally 
invest igatedi^i2i. We flnd that the inelastic part of resonance fluorescence, the part arising 
from nonlinear quantum processes, increases by more than two orders of magnitude with 
respect to the free-dot value. We also investigate the statistics of the photons scattered 
by this system. We show that a small variation of the excitation frequency determines a 
variation of the second-order correlation function for scattered light beyond three orders 
of magnitude. These striking results demonstrate that the arrival of one photon is able to 
suppress or greatly enhance the scattering properties of the second one. 

II. QUANTUM DESCRIPTION OF THE METALLIC NANOPARTICLE- 
QUANTUM DOT INTERACTION 

We consider a spherical QD interacting with a spherical MNP of radius separated by 
a distance R. The system setup is described in Fig. la. The entire system is embedded in 
a dielectric medium with constant permittivity eb- The QD dot is described as a two level 
quantum emitter with a dipole moment /x. This artiflcial molecule is excited by an applied 
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electromagnetic field Ei = £"06"^^^^ + c.c. polarized along the system axis. The positive 
frequency component of the electric field oscillating as exp(— icjt) felt by a quantum emitter, 
E'qd = Eq + is due to the superposition of the free field input field Eq and the field 
arising from the induced polarization of the MNP—: 

where = 2 for an applied electrical field parallel to R {sa = —1 for a field orthogonal to 
R). Pm is the polarization of the MNP, 

P„^4.Wr3(s. + jJ^£|^). (2) 

Fx = is the positive frequency component of the QD polarization, being (a) the ex- 
pectation value of the lowering transition operator a = \g){e\. The frequency dependent 
complex coefficient /3 = {e^{uj) — eb)/(2eb + determines the SP dipole resonance fre- 

quency cjgp satisfying Re[6m(cc;sp)] = —26b. Performing the first order expansion of Re[em(cL;)] 
around cjgp, /3 can be well approximated by the complex Lorenzian 

(3) 



with 

fd R. 

n = I 

duj 

and 7sp = 2r]lm[e^{(jJsp)]. This approximation enables the description of SP resonances 
within the quasi-mode approach, largely exploited in the framework of cavity QED^*^. 
Within this scheme, it is possible to describe the SP field as single resonant mode with a 
linewidth arising from the interaction with a reservoir—. The full quantum dynamics of 
the coupled nanosystem can be derived from the following master equation for the density 
operator, 

p=^[p,Hs]+C,^ + C^. (5) 
The system Hamiltonian is Hs = Hq + H^^t + ^drive with 

Hq = hujsp a^a + h(jj^ a^a , (6) 

being a the Bosonic destruction operator describing the SP field mode and huj^ the energy 
of the QD excitonic transition. The Hamiltonian term describing the interaction between 



the QD exciton and the quantized SP field, in the rotating wave approximation reads 

^^int = ifig {d^cr - aa^) , (7) 

where Tig = jiE^ being iSa = the positive- frequency electric field operator at the QD 
position. The system excitation by a classical input field is described by the following 
Hamiltonian term, 

^drive = -^0 {Xa^ + - M^O {(T^ + ^) • (8) 

The Markovian interaction with reservoirs determining the decay rates 7x and 7sp for the 
QD exciton and the SP mode respectively, is described by the following Liouvillian terms^S, 

^sp = -^{2apa^ — a^ap — pa^a) , (9) 

L^ = — {2a pa^ - a^ap - paV) . (10) 

Starting from [5l the coupled equations of motion for the SP-field expectation values (a) = 
Tr[ap] and for the emitter transition-operator determining the QD polarization (a) = Tr[(jp] 
can be obtained, 



«(wsp-cc;) + ^ («) + + ^^Y^ , (11) 



It 



o) = -\ (00. - ^) + y] (^) - M + 2^ WV) + ^ (1 - 2{a^a)) . (12) 

The equation of motion [12] for the exciton operator expectation value is coupled to higher 
order expectation values. It is known in cavity QED that master equations as ([5j) induce an 
open hierarchy of dynamical equations which needs some approximation. A widely adopted 
truncation scheme is the one based on the smallness of the excitation density which allows to 
truncate with respect to the number of photon number-states to be included^^. The inclusion 
of one-photon states only gives linear optical Bose-like dynamics. Nonlinear optical eflFects 
need the inclusion of multi-photon states. This truncation scheme provides accurate results 
also for strongly interacting systems provided that a suflScient number of photon-states are 
included. At steady state [H] can formally be solved as. 



Equating the obtained electric field expectation value = i£{a) with the corresponding 
field in [H we obtain, 



(14) 



and 



X = el,^/l27^hr]eor^. (15) 
Figure lb displays the dependence of g on the QD-MNP distance R. Throughout the paper 
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FIG. 1: (a) Interaction between a quantum dot and a silver nanosphere: the applied electromagnetic 
field induces a polarization that causes dipole-dipole coupling. States \g) and |e) are coupled via 
the localized surface plasmon dipole mode with a strength g. (b) Dependence of the coupling g on 
the metallic nanoparticle-quantum dot distance . (c) Contour plot of the effective damping rate 
of the quantum dot due to the presence of the metallic nanoparticle as a function of distance and 
detuning. 

we use a dipole moment /x = ero with tq = 0.7 nm (corresponding to 33.62 Debye), being 
e the electron charge. Throughout the paper we use 6b = 3 and we consider a silver MNP 
whose frequency-dependent dielectric permittivity is taken from Ref.^''. Inserting [13] into 
the second term on the r.h.s. of [12] , we derive the effective damping rate of the QD due to 
the presence of the MNP: 
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(16) 



Fig. Ic shows the dependence of Fqj^ on the QD-MNP distance and energy detuning. FinaUy 
comparing the expression for the MNP polarization [2] with [131 we obtain 

= X{0). (17) 

After the determination of £ and X) El estabhshes a precise theoretical framework for non- 
perturbative SP-QED. In this way the dynamics of the emitter interacting with the field 
mode can be treated beyond the Wigner-Weisskopf approximation^-. 



III. THE MNP-QD MOLECULE AS A SATURABLE SCATTERER 

We calculate the Rayleigh scattering adopting the standard method which is valid when 
the size of the scattering objects is much smaller than the wavelength of incident lighlp^. 
In this case the scattering intensity is proportional to Is = {P~P^)^ where = + M^r 
is the total polarization operator and P~ = {P^Y . It is worth noticing that the scattered 

2 

intensity contains a coherent part 1^^^ = (P^) as well as incoherent contributions with 
the frequency of the scattered photons not necessarily coincident with that of incident light 
jmcoh ^ = Ig — /^^^. Figure 2a displays scattering spectra as function of the 

frequency of the incidence light obtained for diflFerent QD-MNP distances R as indicated in 
the panel. The spectra in Fig. 2a have been calculated in the limit of very low excitation 
intensity, where the excitonic populations (crV) <C 1. At i? = 14 nm a Fano-like lineshape 
around the QD transition energy (jJq is evident. For a particular input frequency the scattered 
light is higly suppressed, while at slightly lower energy an enhancement of scattering due to 
constructive interference can be observed. For comparison the plot at i? = 14 nm shows the 
scattering spectrum in the absence of the QD (dash-dotted line). Increasing the distance 
i?, the Fano resonance narrows, due to the reduction of Fqj^. While at i? = 18 nm the 
destructive interference remains almost complete, at larger distances {R = 25 nm), the 
Fano interference effect lowers and the suppression as well as the increase of the scattered 
light are reduced. Figure 2b displays scattering spectra obtained for QDs with different 
excitonic energy levels. Interestingly, the interference effect determining a strong suppression 
of scattering at specific wavelengths of the input field requires no special tuning unlike 
analogous effects in cavity QED^^. While SPs supported by the MNP can be described as 
harmonic oscillators, the single QD displays nonlinearities at single photon level. Hence we 
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FIG. 2: (a) (Color online) Scattered light intensity spectra (red continuous line) calculated for 
different QD-MNP distances R at low density excitation power. For i? = 14 nm, the scattered 
light without the presence of the QD (dot-dashed line) is also plotted, (b) Scattered light intensity 
spectra calculated at i? = 14 nm. Each panel shows calculations for a specific exciton-SP energy 
detuning indicated in the figure. The black continuous line describes plots obtained for an input 
intensity field = 0.02 meV, the red short-dashed line plots obtained at = 0.4 meV, and the 
blue short-dotted line plots at = 1 meV. 



may expect important nonlinear optical effects when increasing the excitation power. Panel 
2b puts forward the dependence of light scattering on the intensity of the input field. The 
continuous lines describe low-field spectra obtained for a Rabi energy ft = 2/x£'o = 2 x 10~^ 
meV. Increasing the input field to Q = 0.4 meV, saturation effects appear (dashed line). At 
Q = 1 meV saturation is almost complete. The hybrid artificial molecule thus behaves as a 
frequency dependent saturable scatterer. 
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IV. RESONANCE FLUORESCENCE 



As outlined in the previous section, when hght excites resonatly a quantum emitter, the 
output scattered hght may be separated into coherent and incoherent parts. The coherent 
part is due to the elastic Rayleigh scattering and has exactly the same wavelength of the 
input light: (jJ^ = cdi- The incoherent part origins from the combination of nonlinear in- 
teractions and quantum fluctuations and cannot be described in a semiclassical context . 
For low incident light intensities elastic scattering dominates, whereas for higher incident 
light intensities two photon processes with + uj^f = 2(jJ\ take place with the energies of 
scattered photons uo^ and uo^' not necessarily equal to uj\. We are interested in calculating 
the steady-state fluorescence emission spectra: 

poo 

Siuj,) = lim 2Re / (p-(t),P^(t + r)) e'^^^dr . (18) 

According to the quantum regression theorem^, two-time correlations {An{t)Ajn{t + r)) fol- 
low the same dynamics of one-body correlation functions (^^(r)) but with the one-time cor- 
relation {An{t)Ajn{t)) as initial conditions. The controlled enhancement of single-molecule 
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FIG. 3: (Color online) (a) Incoherent intensity emission as function of the driving field intensity 
obtained for different exciton-SP detunings. The black short-dotted hne describes the QD inco- 
herent emission in the absence of the MNP. (b) Resonance fluorescence spectra: intensity of the 
scattered field as a function of the dected frequency calculated at three different pump intensities. 

fluorescence due to its near-fleld coupling with a metal nanoparticle, acting as a nanoan- 
tenna, has been widely demonstrated. In this section we investigate the strong influence of 
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such a nanoantenna on the nonhnear quantum scattering of an individual quantum emit- 
ter. Figure 3a displays the incoherent part of the emitted intensity as a function of the 
input field for three different excitonic energy levels. Calculations have been performed for 
a frequency of the incident field uji = uj^. The obtained results show that the MNP is able 
to strongly infiuence resonant quantum optical nonlinear processes. Calculations have been 
performed considering a QD-MNP distance i? = 14 nm. At resonance {uj^ = cjgp) we find an 
enhancement of a factor ^ 260 with respect to the QD incoherent emission in the absence 
of the MNP. The enhancement is the result of two different competing processes: (i) the 
enhancement of the QD emission due to the SP-increased density of photon modes; (ii) the 
quenching of the emitted light due to the losses induced by the MNP described by the rate 
Fqj^ . Fig. 3a shows that, despite the MNP enhancement of the driving field felt by the QD, 
its presence increases the threshold Rabi energy for the process. This is a direct consequence 
of the losses increase induced by the coupling of the QD with the MNP. Figure 3(b) dis- 
plays the resonance fiuorescence spectra (the incoherent part) as function of the detection 
frequency for different excitation powers. It is seen that, with the increasing driving field 
intensity, the single-peak spectrum around u) rm = is transformed into the three-peak 
MoUow spectrum. The emission spectra, owing to the MNP induced losses are strongly 
broadened with respect to that of the free QD. This should enable a better filtering out of 
the strong elastic component which is necessary in order to measure the MoUow spectrum. 



V. PHOTON CORRELATIONS 



The strongly nonlinear atomic response at the single-photon level leads to pronounced 
modification of photon statistics that cannot be captured by only considering average intensi- 
ties, but appears in higher-order correlations of the emitted and scattered fields. Specifically, 
we focus on the normalized second-order correlation functions g^'^'^ for the scattered field 
which for a stationary process can be expressed in terms of the total polarization operators 
as 

g^^\r) = {P-i^P-ii + r)PHt + r)Pnt)) ^^^^ 

Figure 4b shows the normalized second-order correlation functions for scattered photons 
calculated at low excitation power for two specific frequencies of the driving field indicated 
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FIG. 4: (Color online) (a) Comparison between two photon detection probability (short-dashed red 
line) and the squared scattered light intensity (continuous black line) calculated at low excitation 
power, (b) Normalized second-order correlation functions for the scattered field calculated at 
low excitation power for two specific frequencies of the driving field indicated by arrows in the 
inset. The inset displays a detail of the low-excitation power scattering spectrum obtained for 
A = cjx — ^sp — —60 meV. (c) Normalized second order correlation function for the incoherent 
part of the emitted field, at high (black continuous line) and low (red short-dotted line) excitation 
density, (d) As (b) but at higher excitation density. 

by arrows in the inset. The inset displays a detail of the low-excitation power scattering 
spectrum obtained for cjgp — c^x = 60 meV. The continuous line in Fig. 4b shows a huge 
bunching effect. It describes the behavior of an efficient single-photon switch. Scattering 
of single photons at this frequency is highly suppresses due to destructive interference, but 
the arrival of the first photon saturates the QD transition thus enabling the scattering of 
the second one. For delay times among the two detection events larger than the exciton 
decay rate, g^'^^r) 1, which is the standard level for coherent classical light. The one 
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that would be measured in the absence of the QD. The dashed hne in Fig. 4b describes the 
second-order correlation function g^'^^r) obtained fixing the frequency of the driving field at 
the Fano-peak corresponding to constructive interference. In this case g^'^\r) drops below 
the classical level and displays an antibunching effect. Second order correlation functions 
of the electromagnetic field below 1 cannot be describe by classically and are a signature of 
the quantum nature of light. The MNP-QD system thus is able to affect dramatically the 
photon-statistics of scattered light. A small variation of the excitation frequency determines 
a variation of the second-order correlation function for scattered light beyond three orders of 
magnitude. This striking effect can be better understood comparing the expectation values 
for two detection events (the numerator of g^'^\0)) and the squared scattering intensity (the 
denominator) separately as shown in Fig. 4a. {P~ (t) P~ (t) (t) (t)) displays a reduction 



of interference effects with respect to 



. Only one of the two detected photons 



(p-(t)p+(t)) 

is affected by interference: the first photon saturates the QD transition, causing the absence 
of any interference effect on the second one. The second order correlation function for higher 
excitation power is displayed in Fig. 4d. Increasing the excitation power saturates the QD 
and photon statistics approach that in the absence of the QD. Fig. 4c displays the second 
order correlation function only for the incoherent part of the emitted field. This field origins 
from the QD only and its second order correlation function can be calculated as follows, 

Oil^ir) = {aHt)a\t + r)a{t + r)a{t))/\{aHt)a{t))\' . (20) 

9'incoh(0) drops to zero as in standard resonance fiuorescence in the absence of the MNP. After 
the emission of the first photon, the quantum emitter is in its ground state and it needs 
time in order to get excited once again. Interestingly the presence of the MNP determines 
a strong decrease of the excitation time. After r = 2 ps 5'incoh('^) ^ ^- -^^S- 4c has been 
obtained for uji = (jU^ = ^sp- This strong decrease of the excitation time can be very useful 
since it reduces the switching time allowing for ultrafast quantum processing. 



VI. CONCLUSIONS 



We studied theoretically the quantum optical properties of hybrid molecules composed 
of an individual quantum dot and a metallic nanoparticle. The presented numerical calcu- 
lations have been carried out including quantum fiuct nations at all orders. The infiuence of 
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the MNP on the QD has been treated nonperturbatively beyond the Weisskopf-Wigner ap- 
proximation. We apphed the developed theoretical framework to investigate the resonance 
fluorescence of this hybrid system. Its incoherent part, the one arising from nonlinear quan- 
tum processes, results to be enhanced by more than two order of magnitude as compared to 
that in the absence of the metallic nanoparticle. Scattering spectra at different excitation 
powers and calculations of intensity-field correlation functions show that this system can act 
as a nonlinear two-photon switch for incident photons, where the presence (or absence) of 
a single incident photon field is suflicient to allow (or prevent) the scattering of subsequent 
photons. Our calculations have also shown that a small frequency shift of the incident light 
fleld may cause changes in the intensity field correlation function of orders of magnitude. 
The ultracompact hybrid system here investigated is in experimental reach. The intriguing 
quantum optical nonlinear properties here demonstrated provide indications that these hy- 
brid systems could be used as ultra-compact elements in quantum-information technology, 
and for single-photon devices. The nonperturbative method here developed to study the 
coupling of an individual optical emitter to the quantized SP field, can be extended to study 
the quantum optical properties of more complex structures as a trimer composed by two 
MNPs and a QD placed in their gap, displaying the vacuum Rabi splitting—. 
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